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Abstract 

The Schwinger-boson mean-field theory is used to study the three- 
dimensional antiferromagnetic ordering and excitations in compounds 
L2BaNiO^, a large family of quasi-one-dimensional mixed-spin antiferromag- 
net. To investigate magnetic properties of these compounds, we introduce a 
three-dimensional mixed-spin antiferromagnetic Heisenberg model based on 
experimental results for the crystal structure of L2BaNiOrj. This model can 
explain the experimental discovery of coexistence of Haldane gap and antifer- 
romagnetic long-range order below Neel temperature. Properties such as the 
low-lying excitations, magnetizations of Ni and rare-earth ions, Neel temper- 
atures of different compounds, and the behavior of Haldane gap below the 
Neel temperature are investigated within this model, and the results are in 
good agreement with neutron scattering experiments. 

PACS numbers: 75.10.Jm,75.50.Ee,75.50.-y 
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I. INTRODUCTION 



The existence of Haldane gap [Q] in the magnetic excitation spectrum has kept the integer- 
spin one-dimensional (ID) Heisenberg antiferromagnet to be one of the most interesting sub- 
jects in condensed matter physics during the past twenty years. The integer-spin Heisenberg 
antiferromagnetic (AF) chain should have a singlet ground state, exponentially decaying cor- 
relations, and a quantum gap, which have been confirmed by numerous theoretical [^^^,^,|^ 
and experimental 0,^,^|TU| studies. Kennedy and Tasaki have proved that the appearance 
of Haldane gap in the spin-1 AF chain corresponds to the breaking of a hidden x Z2 
symmetry fill]. The recent discovery of coexistence of Haldane gap and AF long-range order 
(AF LRO) in rare-earth compounds L2BaNiO^ with L=Y, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, 
and Tm |[T0|JT^JT3|JT^JT5|JT^ has offered an opportunity to investigate the effect of staggered 
field on the Haldane chain. The polarized and unpolarized inelastic neutron scattering ex- 



periments on L2BaNiO^ |ll^ , p!6[| present that, on the one hand, there are AF interactions 
between spins of rare-earth and Ni ions, which lead to the three-dimensional (3D) AF LRO 



with Neel temperature {Tn) ranging from 20 — 70/^ |T^|T^. On the other hand, these com- 
pounds are also characterized by the presence of NiO^ octahedron chains along a-axis. As 
is anticipated, above the Neel temperature, energy gap has been found in the excitation 
propagating along A^i-chain. Also the Haldane gap is discovered to persist in the 3D AF 



phase and increase with decreasing temperature below T/v [|14|,|T5|,|16 



Recently several theoretical works have focused on the quasi-lD Haldane system 
L2BaNiO^ [ll7| , [T8| , [T9| . To describe the coexistence of Haldane gap and AF LRO, the non- 
linear sigma model in a static staggered magnetic field has been introduced by Maslov and 



Zheludev [|T^. This system has been studied numerically within the density matrix renor- 
malization group method by adopting the model of an ID spin-1 Heisenberg chain in a static 



field IjT^ . In the above studies, the reason for introducing the static staggered field is based 
on the assumption that there are 3D directly AF exchange interactions between the spins 
of rare-earth ions, which lead to the AF LRO below Tjq. However, it has been pointed out 
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by some experiments that spins of Ni ions also play an important role in forming 

AF LRO, and the two-dimensional mixed-spin Heisenberg model with Si = 1 and S2 = 1/2 
chains staked alternatively has also been adopted for this system [|T^. 

Depending on the experimental finding of crystal structure, we introduce a more suitable 
3D mixed-spin AF Heisenberg model to describe the magnetic properties of compounds 
L2BaNiO^. In this model we define si = 1 for the spins of Ni ions, and investigate different 
members of this family by changing the spin value S2 of rare-earth ions correspondingly. 
We use the Schwinger-boson mean- field (SBMF) theory to study this 3D mixed- 

spin model. The SBMF theory has been successful to the study of spin-integer Heisenberg 
chains p0|j21[| , and has also been extended to cases with magnetic ordering by identifying 



the magnetization with the Bose condensation of the Schwinger bosons [^,^,^. Within 
this model, we could explain the experimental discovery of coexistence of Haldane gap and 
AF LRO. The interactions between spins of Ni and rare-earth ions are proved to be very 
important in forming the AF LRO. Properties such as the magnetization, Haldane gap and 
Neel temperatures are also discussed for cases with different 82- Comparing our results with 
neutron scattering experiments, we find that our findings could explain some experimental 
discoveries. 

The paper is organized as follows: In Sec. II we introduce the 3D mixed-spin Heisenberg 
model and Schwinger-boson mean-field theory. Our results of the magnetic properties of 
L2BaNiO^ compounds are presented in Sec. HI. Finally, we conclude our findings in Sec. 
IV. 

II. MIXED-SPIN HEISENBERG MODEL AND SCHWINGER-BOSON 

MEAN-FIELD THEORY 

The crystal structures of L2BaNiO^ compounds have been investigated extensively by 



neutron scattering experiments |12,13,r^. They belong to the orthorhombic system, having 



approximate cell Parameters around a = 3.8A, h = 5.8A and c = 11. 3A. In all members 
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of this compounds, as figure 1(b) shown, strings of distorted NiO^ octahedron share apical 
oxygen and form Haldane chains along the a-axis. The intrachain Ni-O-Ni AF superex- 
change coupling is about 200-300A' |T2|JT3| , p!^ , which is the strongest magnetic interaction 



in the compound. While in the he plane, as figure 1(a) shown, there is one oxygen between 
the Ni and rare-earth (L)ions, and thus the Ni-O-L interactions establish links between 
individual Ni chains. This coupling is of the order of tens of Kelvin (lO-SOfC). Figure 1(a) 
also shows that there are L-O-L AF interactions between rare-earth ions, which is of the 
order of IK or less. 

Based on the experimental discovery of the crystal structure of compounds L2BaNiO^, 
we introduce a 3D mixed-spin AF Heisenberg model expressed by the Hamiltonian 

i,Va i,Vb,Vc j 

where 5*1 and S2 are the spin operators of Ni and rare-earth ions respectively, rjb and rjc 
denote a sum over nearest-neighbor (NN) bonds in the be plane, and so is rja for the NN 
bonds along the a-axis. The AF superexchange couplings along the Ni chain, between L 
and Ni ions, and among the rare-earth L ions are represented as Ji, J2, and J3 respectively 
as figure 1(a) shown. 

The Schwinger-boson theory introduces Si = ^bf^aa^fshis (a (or /5)=t, i) [^,^, and the 



spin degrees of freedom are mapped to the boson degrees of freedom. Meanwhile the original 
spin Hilbert space corresponds to a boson Hilbert subspace in which bf^bi-^ + = 2sj. 
These constraints to the boson Hilbert space are imposed in the Hamiltonian (1) by introduc- 
ing two kinds of Lagrangian multipliers Xi{i) and A2(j). In addition, we define the bond op- 
erators as Qa{i,Va) = felit&l(i+r,„)i-fclii&l(i+r,„)T, Qhii,Vb,Vc) = 6liT&2(i+r,6+,,,)i " &ia&2(i+r,5+T,,)T , 

and Qc{j, Co) = ^2jT^20+c?))i ~ ^2ii&2(i+c^))f "^^^ Hamiltonian (1) can be rewritten as 

H=-IJiY: {Qaih Va)Qa{l, Va) " 2s?} - ^ J3 ^ {g+(j, Co)a(j, Cq) - 2sl} 

'i.,Va j 

"^•^2 {Qh{hVb,Vc)Qh{i,Vb,Vc) -'2S1S2} 

i,Vb,Vc 

+ E {Ki^bu^ + bt^M - 2si} + E Hj) {btjM + btjib2ji - 2S2} . (2) 



Next, we make a Hartree-Fock decomposition of Eq. (2) by taking the average values of 
the bond operators and Lagrange multiphers to be uniform and static as {QaihVa}) = Qa, 
{Qhi'i,'nb,Vc)) = Qh, {Qc{j,co)) = Qc, (Ai(i)) = Ai, and (A2(j)) = A2. Under the Fourier 
transformation, we obtain the following mean-field Hamiltonian in the momentum space 

^MF ^Y^2Zk {bik^b^-k)i - bikibi(-k)^ + h.c} + ^4Dk {6ifcT^-fe)i - ^ifei^-fc)T + h.c} 

k k 

+ Y^'^Xk {&2fcT^2(-fc)i - &2fei&2(-fe)T + h.c\ + [4:\iblf.^bika + 8A262fea^2fea} 
k ka 

+2N {jiQl + AJ2QI + J3QI - 4AiSi + 2Jisl + 8J2S1S2 + 2 Jgs^ - 8A2S2} , 

where Zk = JiQaCOs{k^ao), Xk = JsQcCOsik^Co), and Dk = 2J2QhCOs{k^Co) cos{kybo). 
Diagonahzing in Eq. (3) by the Bogoliubov transformation, we obtain 

//^^ = 2 ^ {E+ {atk,a,ka + OC^kMka) + {PtkaPlka + P2kaPtka) } + ^0 (4) 
k 

with 



(3) 



lA±VA^-4:B 



Ei-- ^ 2 
A = Xi + A\l^ [Zl + xl + '^Dl^ 
B = 4A1A2 (A1A2 - Dt) - xl^l - ^ZlXl + (xkZk - 
Eo = 2N {jiQl + A.hQl + hQl + 2.hs\ + 8J2S1S2 + 2J3S2} 

-SiVAi f si + ^"j - 16iVA2 {s2 + ^) . (5) 



2J V 2 

Furthermore, the mean-field free energy is given by 

= I X: In { (2 sinh ) + In (2 sinh ) | + (6) 

The mean-field equations are obtained by differentiating F^^^^ with respect to the pa- 
rameters Qa-i Qh, Qc, Ai, and A2 as 

QpMF QpMF gpMF QpMF QpMF 



= 0. (7) 



dQa dQh dQc dXi 8X2 
Thus we obtain five self-consistent equations to determine the average values of the bond 
operators and Lagrange multipliers. 



III. MEAN-FIELD SOLUTIONS 



In this section, we present the solutions of the SBMF theory. The experimental discovery 
of coexistence of Haldane gap and AF LRO below T/v are obtained by our study. As equation 
(5) shown, there are two branches of the magnetic excitations and , which have quite 
different behaviors below T^. On the one hand, the magnetic excitation E^ has an energy 
gap in the whole temperature region, and below T/v this energy gap increases with decreasing 
temperature. On the other hand, within the temperature region from zero to T/v, E^ keeps to 
be gapless and has its minimal value E'^ = at A; = 0. Under this condition, the Schwinger- 
boson condensation occurs and leads to the AF LRO in this system. In our calculation, we 
introduce the Schwinger-boson condensation into the self-consistent mean-field equations, 
and obtain the temperature dependence of the magnetizations of Ni (Mi) and rare-earth 
(M2) ions below T/v respectively. The magnetizations Mi and M2 are expressed as 



2N 

M2 = (5f ) = ^ - btn^m)- (8) 

j 

We choose the AF superexchange interactions as Ji = J, J2 = O.IJ and J3 = O.OIJ (J > 0) 
according to the experimental studies of the magnetic properties of compounds L2BaNiO^ 



T^,[T3|jTj]. The AF interactions between Ni ions has been estimated as 2 J = 200 — 300fC by 



the neutron scattering experiments []12|,|14[. To make it simple and clear, in our calculation 
we choose J = lOOi^. 

The neutron scattering experimental results of Nd2BaNiO^ obtained by Yokoo et al. 



T5| are shown in Fig. 2(a) and 2(b) (open circles). The solid lines in Fig. 2(a) and 2(b) 
show fits to the above experimental results by our numerical calculations for the case with 
S2 = 3/2 respectively. The temperature dependence of magnetizations Mi/si and M2/S2 
for the cases of S2 = 1/2, 1, 3/2, 2, 5/2, and 3 are shown in figure 2(c) and 2(d) respectively. 
As temperature increases, the thermal fiuctuation in the system becomes stronger, and the 
magnetizations decrease rapidly and drop to zero at the Neel temperatures. Therefore, we 
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could also determine the Neel temperatures Tn for all members of the L2BaNiO^ family 
through our calculation. 

In figure 3, we plot the Neel temperature T/v as a function of spin value S2 (filled circles 
and dotted line). For comparison, the Neel temperatures of compounds Ho2BaNi05 (T/v = 
53K) ig, Nd2BaNiO^ {Tn = and Pr2BaNiO^ {Tn = 2AK) obtained 

by the neutron scattering experiments are also shown in Fig. 3 (filled triangles). Our 
theoretical results roughly agree with the experimental results. We also obtain that the 
Neel temperature Tn increases monotonously with the increasing of S2- As the spin value 
of the rare-earth ions S2 equal to 1/2 and 3, we get that the minimize and maximum Neel 
temperature are T^^*" 0.207J and T^'"-^ 0.872J respectively. Our result of the Neel 
temperature region is approximately from 20. 7/^ to 87.2K when J = lOOK, which is in 



good agreement with the experimental estimation of Tn region 20 — 70K [[T2| , p!3[ . 

The two branches of the magnetic excitation and Ej^ represent the spin fiuctuations 
along the a-axis and within be plane respectively. In the exactly one-dimensional case ( J2 = 
J3 = 0), the excitation Ej^ vanishes and the energy gap of E'^ is just the Haldane gap of 
AF Heisenberg chain, which is closely related to the breaking of a hidden Z2 x Z2 symmetry 



|rT|. We obtain that the coexistence of Haldane gap and AF LRO below Tn is a common 



feather of all members of compounds L2BaNiO^ except Y2BaNiO^. The neutron scattering 
experiments have discovered that, below the Neel temperature, the Haldane gap increases 
as the temperature decreases |]14|,|15|,|16| . The temperature dependences of the Haldane gap 
A for the cases of S2 = 1/2, 1,3/2,2,5/2 and 3 are also investigated by the SBMF theory 
and the behaviors are shown in figure 4(a) respectively. The temperature dependence of the 
energy gap A obtained by our calculation is found to agrees with the experimental discovery. 

Our calculation also implies that the effect of the staggered magnetization on the mag- 
netic excitation is to widen the Haldane gap. Below T/v, the effective internal magnetic 
field Heff imposed on the Haldane chain is assumed approximately as the magnetization 
Ml of Ni ions. This field is found to increase with the decreasing of temperature and the 
increasing of spin value S2 because of the thermal fiuctuation being weakened and the AF 



ordering being enhanced. In figure 4(b), we plot the energy gap A of zero temperature as a 
function of H^ff, and we find that -f^e// has strong effect to widen the Haldane gap. 

In addition, based on the 3D mixed-spin model, we obtain that AF LRO below the 
Neel temperature is not constructed only by the rare-earth ions. Our results support the 
suggestion that Ni ions also play an important role in forming the AF LRO. We plot in figure 
5 the Neel temperature as a function of AF coupling J2 in the cases with S2 = 1/2, Ji = J, 
and J3 = O.OIJ. In compounds L2BaNiO^, the effective interactions between individual Ni 
chains rely on the AF coupling J2 between Ni'^'^ and rare-earth L^"*". We obtain that the 
Neel temperature rises rapidly with the increase of J2 as shown in figure 5, so the coupling 
J2 are important in forming the AF LRO. Besides, as J2 = 0, there are no interactions 
between Ni chains and thus the Neel temperature drops to zero. The Haldane excitation 
energy in 3D and ID cases are shown in figure 6(a) and 6(b) respectively. Here we choose 
S2 = 1/2, Ji — J, J2 — O.OIJ, J3 = O.OIJ, and kx — ky — for the 3D case, and obtain that 
the corresponding Neel temperature is — 0.206 J. To compare the behaviors of magnetic 
excitations below and above the Neel temperature, we study two conditions of T = 0.1 (solid 
hues) and T — 0.4J (dotted hues) respectively for both ID and 3D cases. In the 3D case, 
we find, in figure 6(a), that the energy gap below the Neel temperature is obviously bigger 
than that above the Neel temperature, which is in opposition with the behavior for pure ID 
case (shown in figure 6(b)). In addition, we obtain that the thermal fluctuation has strong 
effect to destroy the 3D spin correlations in the compounds, as a result the behaviors in 3D 
case is the same as that in ID case when temperature is above the Neel temperature. 



IV. SUMMARY 

In conclusion, we have introduced a 3D mixed-spin AF Heisenberg model based on the 
experimental results of the crystal magnetic structure of compounds L2BaNiO^, and studied 
this model with the SBMF theory. The experimental discovery of coexistence of Haldane 
gap and AF LRO below TV has been deduced by our calculation. Properties such as the low- 
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lying excitations, magnetizations of Ni and rare-earth ions, Neel temperatures of different 
members of this family and behavior of Haldane gap below TV have also been investigated 
within this model. 

We have obtained two branches of the magnetic excitations E'^ and , of which E^ 
has an energy gap and E^ is gapless below Tjq. The theoretical result of the Neel temper- 
ature region is approximately from 20. to 87.2K, which is in good agreement with the 
experimental estimation of the region 20 — 70K. We have also obtained that Haldane gap 
increases with decreasing temperature, and the effect of the magnetization is to widen the 
Haldane gap. Our results are in good agreement with the experimental discoveries. Our 
findings also support the suggestion that the AF LRO below Tjv is not constructed only by 
the rare-earth ions, and Ni ions also play an important role in forming the AF LRO. 

ACKNOWLEDGMENTS 

We thank Dr. Xintian Wu for helpful discussions. This work was supported by the 
National Natural Science Foundation of China under Grant Nos. 10125415, 10074007, and 
90103024. 



9 



REFERENCES 

[1] F. D. M. Haldane, Phys. Lett. 93A, 464 (1983); Phys. Rev. Lett. 50 , 1153 (1983). 

[2] I. Affleck, J. Phys. Condens. Matter 1, 3047 (1989). 

[3] M. P. Nightingale and H. W. J. Blote, Phys. Rev. B 33, 619 (1986). 

[4] K. Nomura, Phys. Rev. B 40, 2421 (1989). 

[5] T. Sakai and M. Takahashi, Phys. Rev. B 42, 1090 (1990). 

[6] S. V. Meshkov, Phys. Rev. B 48, 6167 (1993). 

[7] R. M. Morra, W. J. L. Buyers, R. L. Armstrong, and K. Hirakawa, Phys. Rev. B 38, 543 
(1988). 

[8] S. Ma, C. Broholn, D. H. Reich, B. J. Sternheb, and R. W. Erwin, Phys. Rev. Lett. 69, 
3571 (1992). 

[9] A. Zheludev, S. E. Nagler, S. M. Shaprio, L. K. Chou, D. R. Talham, and M. W. Meisel, 
Phys. Rev. B 53, 15004 (1996). 

[10] D. J. Buttrey, J. D. SuUivan, and A. L. Rheingold, J. Sohd State Chem. 88, 291 (1990). 

[11] T. Kennedy and H. Tasaki, Phys. Rev. B 45, 304 (1992); Rev. Math. Phys. 6, 887 (1994). 

[12] E. Garcia-Matres, J. Rodriguez-Carvajal, J. L. Martinez, A. Salinas-Sachez, and R. Saez- 
Puche, Sohd. State. Commun. 85, 553 (1993). 

[13] V. Sachan, D. J. Buttrey, J. M. Tranquada and G. Shirane, Phys. Rev. B 49, 9658 (1994). 

[14] A. Zheludev, J. M. Tranquada, and T. Vogt and D. J. Buttrey, Phys. Rev. B 54, 6437 
(1996). 

[15] T. Yokoo, et ai, Phys. Rev. B 58, 14424 (1998). 

[16] S. Raymond, T. Yokoo, A. Zheludev, S. E. Nagler, A. Wildes, and J. Akimitsu, Phys. 

10 



Rev. Lett. 82, 2382 (1999). 

[17] S. Maslov and A. Zheludev, Phys. Rev. B 57, 68 (1998); Phys. Rev. Lett. 80, 5786 (1998) 

[18] J. Lou, X. Dai, S. Qin, Z. Su, and L. Yu, Phys. Rev. B 60, 52 (1999). 

[19] Y. Takushima, A. Koga, and N. Kawakami, Phys. Rev. B 61, 15189 (2000). 

[20] D. P. Arovas and A. Auerbach, Phys. Rev. B 38, 316 (1988); D. Yoshioka, J. Phys. Soc. 
Jpn. 58, 32 (1989). 

[21] A. Auerbach, Interacting Electrons and Quantum magnetism, ( Springer- Verlag, 1994). 

[22] S. Sarker, C. Jayaprakash, H. R. Krishnamurthy and M. Ma, Phys. Rev. B 40, 5028 
(1989); C. L. Kane, P. A. Lee, T. K. Ng, B. Chakraborty, and N. Read, Phys. Rev. B 41, 
2653 (1990). 



11 



FIGURES 

FIG. 1. (a) Structural relation between Ni and rare-earth sites in the be plane, and (b) Ni and 

apical oxygen form the Haldane chain along the a-axis. 

FIG. 2. The temperature dependence of the ordered moment on the Ni (a) and Nd (b) sites in 
Nd2BaNiO^ sample taken from Ref. 14 (open circles), and the solid lines show fits to the exper- 
imental discovery by our numerical results of the case with s = 3/2 respectively; Magnetizations 
Mi/si (c) and M2/S2 (d) as functions of temperature below Tn for cases of S2 = 1/2, 1, 3/2, 2, 5/2, 
and 3 respectively. 

FIG. 3. Neel temperatures of the compounds L2BaNiO^ for the cases with different spin value 
S2- The circles represent the theoretical results, and the triangles are experimental findings of 
compounds Nd2BaNiO^, Pr2BaNiO^ and Ho2BaNiO^ taken from Ref. 11-13. 

FIG. 4. (a) The temperature dependence of Haldane gap for cases of S2 = 1/2,1,3/2,2,5/2 
and 3 below the Neel temperatures respectively; (b) The Haldane gap as a function of the effective 
internal magnetic field -ffe//- 

FIG. 5. Neel temperature as a function of AF coupling J2 between Ni'^'^ and ions of case 
S2 = 1/2, Ji = J, and J3 = O.OIJ. 

FIG. 6. The branch of Haldane excitation E'^ as a function of kz in (a) 3D with S2 = 1/2, 
Ji = J, Js = O.OIJ and (b) pure ID cases. The solid lines show results of temperature T = O.IJ, 
and the dotted lines present results of T = 0.4 J. 
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